Strength of the null singularity inside black holes 
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We study analytically the Cauchy horizon singularity inside spherically-symmetric charged black 
holes, coupled to a spherical self-gravitating, minimally-coupled, massless scalar field. We show that 
all causal geodesies terminate at the Cauchy horizon at a null singularity, which is weak according 
to the Tipler classification. The singularity is also deformationally-weak in the sense of Ori. Our 
results are valid at arbitrary points along the null singularity, in particular at late retarded times, 
when non-linear effects are crucial. 
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I. INTRODUCTION 

The issue of spacetime singularities - which are known 
to inevitably occur inside black holes under very plausible 
assumptions |ffl] — is an intriguing puzzle of physics. The 
laws of physics, as we presently understand them (e.g., 
classical General Relativity), are presumably invalid at 
singularities. Instead, some other theories (e.g., quan- 
tum gravity), as yet unknown, are expected to take over 
from General Relativity and control the spacetime struc- 
ture. The General Relativistic predictions are neverthe- 
less of the greatest importance, as they reveal the space- 
time structure under extreme conditions in the strong- 
field regime. Of particular interest is the possibility that 
there are two distinct ways in which General Relativity 
can fail at different types of singularities; For one type 
of singularities the failure is through infinite destructive 
effects on physical objects, whereas for the other type the 
failure is through breakdown of predictability. 

Until recently, the only known generic singularity in 
General relativity was the Belinsky-Khalatnikov-Lifshitz 
(BKL) singularity M. According to the BKL picture, 
spacetime develops a succession of Kasner epochs in 
which the axes of contraction and expansion change di- 
rections chaotically. This succession ends at unbounded 
oscillations at a spacelike singularity, which is unavoid- 
ably destructive for any physical object - a strong sin- 
gularity. In the last several years, however, evidence has 
been accumulating that the BKL singularity is not the 
only type of singularity which may evolve in General Rel- 
ativity from generic initial data. 

The new type of singularity forms at the Cauchy hori- 
zon (CH) of spinning or charged black holes. (For a re- 
cent review see The features of this singularity are 
markedly different from those of the BKL singularity: (i) 
It is null (rather than spacelike), (ii) It is weak (according 
to Tipler's classification m); Specifically, the tidal defor- 
mations which an extended physical object suffers upon 
approaching the singularity are bounded. In the case 
of a spinning black hole, the evidence for the null and 
weak singularity has emerged from analytical perturba- 
tive and non-perturbative analyses. In addition, 



the local existence and genericity of a null and weak sin- 
gularity in solutions of the vacuum Einstein equations 
was demonstrated in || . This was more recently demon- 
strated also in the framework of plane-symmetric space- 
times in ||. For the toy model of a spherical charged 
black hole, the main features of the CH singularity were 
first found analytically for simplified models based on 
null fluids ]lO|-p^|, and later confirmed numerically for a 
model with a self-gravitating scalar field 13 14]. Expres- 



sions for the divergence rate of the blue-shift factors for 
that model, which are valid everywhere along the CH, 
were found analytically in fl5|] . Those expressions are 
exact on the CH as functions of retarded time. How- 
ever, they are only asymptotic expressions as functions 
of advanced time (see below). 

The strength of the null singularity is of crucial im- 
portance for the question of the hypothetical possibility 
of hyperspace travel through the CH of black holes. A 
necessary condition for this possibility to be realized is 
that physical objects would traverse the CH peacefully. 
Because the CH is known to be a curvature singularity, it 
is necessary that the singularity would be weak according 
to the Tipler classification of singularity strengths. For 
the toy model of a spherical charged black hole, which 
we shall study here, the properties of the CH singularity 
which have been found in |y]-[l5| are all consistent with 
the picture of a Tipler weak singularity. However, the 
weakness of the singularity was demonstrated only for 
the simplified Ori model |L2j and at asymptotically early 
times for spinning black holes g], where there are still 
no strong non-linear effects, such as focusing of the null 
generators of the CH, which are crucial at later times. In 
the context of spherical charged black holes and a self- 
gravitating scalar field, several important features of the 
spacetime structure have been found in fully non-linear 
numerical simulations. Specifically, it was shown that for 
any point along the CH singularity there existed coordi- 
nates for which the metric coefficients were finite and 
the metric determinant was non-degenerate in an open 
neighborhood to the past Jl^Jl^]. However, despite pre- 
vious claims ||[L4|,[L6) , this still does not guarantee that 
the singularity is weak in the Tipler sense |17|. 
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It is the purpose of this paper to present an analyt- 
ical demonstration of the weakness of the singularity 
for the model of a spherical charged black hole with a 
self-gravitating, minimally-coupled, massless, real scalar 
field. Our results are valid at arbitrary points along the 
CH singularity, in particular at late times, where strong 
non-linear effects (focusing of the null generators of the 
CH and growth of the blue-shift factors) are crucial. In 
fact, our results are valid everywhere along the CH sin- 
gularity, down to the event where the generators of the 
CH are completely focused, and the singularity becomes 
spacelike and Tipler strong |l8j. We emphasize that al- 
though our discussion here is analytical, we do make as- 
sumptions which are based on results obtained by numer- 
ical simulations. 

II. STRENGTH OF THE SINGUALRITY 

We write the general spherically-symmetric line- 
element in the form 

ds 2 = -f(u,v)dudv + r 2 (u,v)dfl 2 , (1) 

where dfl 2 = dd 2 + sin 2 $ dip 2 is the line-element on the 
unit two-sphere. The coordinates u, v are any outgoing 
and ingoing null coordinates, correspondingly. (Below, 
we shall specialize to a specific choice of gauge, and de- 
fine a particular choice of an ingoing null coordinate.) 
We consider the class of scalar field perturbations which 
is inherent to any gravitational collapse process. These 
are the perturbations which result from the evolution 
of non- vanishing multipole moments during the collapse. 
When these perturbations propagate outwards, they are 
partially reflected off the spacetime curvature and cap- 
tured by the black hole. This process results in a scalar 
field, which at late advanced times decays along the event 
horizon like an inverse power of advanced time. Specifi- 
cally, we assume that the scalar field behaves along the 
event horizon at late times according to $ EH oc (Kv e )~ n 
|^9|-pT|, where n is a positive integer which is related to 
the multipole moment under consideration. (We do not 
consider, however, other possible sources of perturbations 
p2||.) By v e we denote the usual advanced time in the 
Eddington gauge, and r±, k are the outer and inner hori- 
zons and the surface gravity of the latter, respectively, 
for a Reissner-Nordstrom black hole having the same ex- 
ternal parameters as the black hole we consider has at 
late times. We define the dimensionless ingoing Kruskal- 
like coordinate by V = — exp(— nv e ). In the Kruskal 
gauge we denote the metric function g u y by —F/2. For 
this model, it can be shown analytically that at arbitrary 
points along the CH the following relations are satisfied, 
to the leading orders in [— ln(— V)] 1 fL5| : 

r,v = [ In (-V0P 2 - 2 {l + M- M-VT 1 



+ b 2 [- H-V)]- 2 +0[- ln(-V)p 3 } (2) 

*.v = I- ln (-v)r- 1 {i + ci [- H-V)]' 1 

+ c 2 [- ln(-y)p 2 +0[- ln(-y)]- 3 } • (3) 

Here, A = [r + /(2r_)](r + /V_ + r_/r+), and the expan- 
sion coefficients frj and C{ are functions of retarded time 
only. Note that in the limit V — » these are exact ex- 
pressions as functions of retarded time. That is, to the 
leading order in [— ln(— V)] 1 there is implicit depen- 
dence on retarded time through r = r(u), and along the 
CH singularity both ry and $,y are exactly inversely 
proportional to r(u), in the following sense. Consider 
two outgoing null rays, and let one ray be at u — ui, say, 
and the other at u — u 2 . The ratios ry{2)/ry{\) and 
$ j v(2)/$,v(l) approach r(ui)/r(u 2 ) as V — > 0. Tak- 
ing now u\ to be in the asymptotically-early parts of the 
CH, where r(u\) ~ r_, we find that both ry and $y 
are inversely proportional to r(u). As r(u) is monotoni- 
cally decreasing as a function of retarded time along the 
CH, we find that ry and <&y grow monotonically along 
the CH. This growth is a non-linear effect which indi- 
cates the strengthening of the singularity along the CH 
(although the singularity is still weak according to the 
Tipler classification; see below). 

All the non-zero components of the Riemann- 
Christoffel curvature tensor R^pa are given completely 
in terms of the divergent blue-shift factors ry, Qy, and 
the finite quantities r, r„,$ „, and F. Interestingly, 
Rfivpa does not depend on gradients of F . This can be 
understood from the following consideration. The ten- 
sor Rpupa can be written as the sum of the Weyl ten- 
sor, and another tensor which is built from the Ricci 
and the metric tensors (but not involving their deriva- 
tives). In spherical symmetry the Weyl tensor is given 
completely in terms of the mass function, which is de- 
fined by r nT' 1 * = 1 — 2m(u,v)/r + q 2 /r 2 , q being the 
charge of the black hole. In Kruskal-like coordinates the 
mass function m(u, V) = (r/2)(l + 4r. u ry/F) + q 2 /(2r), 
which depends only on r,ry, r u , and F. (The divergence 
of the mass function at the singularity, and consequently 
also the divergence of curvature, is evident from the di- 
vergence of ry and the fmiteness of r, r „ ,and F.) The 
Ricci tensor R^ u — 2$ iM $ jly , and consequently R^ upa is 
independent of gradients of F. 

We find that the components of Rp. vpa which have the 
strongest divergence near the CH are Ry v y v and Ry-dv-a- 
(In the Appendix we list all the non-zero independent 
components of the Riemann-Christoffel tensor.) It can 
be readily shown that 

Rv-dV-o = -r 2 (<&y) 2 = sm^ 2 dR VvVv . (4) 

We denote these two components schematically and col- 
lectively by 1Z. We find that, to the leading orders in V 



2 



and in [— In (— V)] , 

n(v) a v- 2 1 



ln(-V)] 



-2n-2 



(5) 



(The other divergent components of i? M ^ p(T , i.e., 
RuVvV)R#ip#tp,Ru&v#, and R uv v v are proportional to 
the leading order in V to V~ 1 times a logarithmic fac- 
tor.) Because both metric functions r and F have finite 
values at the CH (which is known from the numerical 
simulations of |ll|[l4|] ) , it is easy to show that the depen- 
dence of TZ(V) on V does not change when we transform 
to a parallelly-propagated frame. 

We next find V(t) as a function of affine parameter 
(proper time) r along a general null (timelike) geodesic. 
For general causal geodesies, the geodesic equations are 

iiv = (mr 2 Cl 2 — p^j (6) 

v + [f, v /f)v 2 + 2mr> (fuv + p) /(rf) = (7) 

u + (f, u /f)u 2 + 2mr iV (fuv + p) /(rf) = 0. (8) 

Here, m = 0(1) for radial (non-radial) geodesies, and 
p = 0(— 1) for null (timelike) geodesies. A dot denotes 
differentiation with respect to affine parameter (proper 
time), and tt 2 = $ 2 + sin 2 -dip 2 . The geodesic equations 
can be solved to the leading order in [ln(— V^)] -1 for all 
causal geodesies. This is done by using the field equation 


F.y/F = r, w/r, v + r($,y) 2 /r, v (9) 

to find F.v/F explicitly. Substituting Eqs. (@) and (§) 
in Eq. (^|) we find 

(\nF) v = [In (-r, v )] tV + ^ + (2ci - h) - [- In (-V)]' 1 
+ (b\ - 6 2 - 26ici +c\ + 2c 2 ) i [- In {-V)Y 2 

+ o|l[-in(-y)]- 3 |. (10) 

Integration yields 
lnF = InFo + ln(V>,y) - (2c x - 6 X ) {In [- ln(-F)]} 
+ (6 2 - 6a - 26ic a +c 2 + 2c 2 ) [- In (-F)]" 1 
+ 0{[-ln(-V)}- 2 }. (11) 

Here, lni^o is an integration constant, which can be a 
function of u. Exponentiating both sides, and substitut- 
ing Eq. (0) for r y we find 



From the numerical results of Jl 3|J 1 4| ] it is known that as 
V — > 0, F approaches a finite value. Consequently, in 
order to have consistency with the numerical results we 
require that b\ — 2c± — 2n + 2, which implies that 



F = F( 



F = F 



~ (nr-A) . n -2n-2+6i-2c 1 



ln(-V-)]" 



X {l + h [- In (-V)]" 1 + O [- In (-V)P 2 } 

x exp { (b\ - 6 2 - 26ici + c 2 + 2c 2 ) [- In (-V)f 1 

+ 0[-ln(-V)}- 2 }. (12) 



oH {l + B[- In (-V)]- 1 + O [- In (-V)}- 2 } , (13) 



where B = (2n + 3)6i - 6 2 + c 2 + 2c 2 , and F (u) = 
Fo(nr_A) 2 /r. For the logarithmic derivative of F we 
find, to the leading order in [— In (—V)] , that 



F, v /F = B^[-\n(-V)]- 2 . 



(14) 



(Higher order terms in [— ln(— V)]" are functions of 
retarded time.) We note that only b\ and c\ are 
constrained. The coefficients of highr-order terms in 
[— ln(— V)] 1 are immaterial near the CH for our de- 
termination of the strength of the singularity. Note that 
F — > Fo(u) as V — > 0, and that F is not analytic in 
P. In fact, this is an important property of the CH sin- 
gularity: In a kruskal-like gauge the metric functions r 
and F are finite at the singularity, but their gradients in 
the outgoing direction diverge. The finiteness of r and 
F at the CH also implies that the metric determinant 
is non-degenerate. This expression for F is similar to 
the behavior of the g u y metric function found for the 
simplified Ori model |l2[| . We stress that although this 
expression for F is exactly valid everywhere along the CH 
singularity, it still does not allow us to find the variation 
of F with retarded time, as we do not know the form 
of Fq(u) or r(u) along the CH. We note that near the 
CH the metric function F is monotonic in V . This result 
is in accord with the numerical results of fll4|| , (Notice, 
however, the disagreement with the numerical results of 
Jl6| . It is reasonable to expect the behavior of the metric 
functions near the CH to be similar for both cases of real 
and complex scalar fields. The lower panels of Figs. (3) 
of Jl6[ imply, however, a non- monotonic behavior of F. 
That kind of behavior can be obtained from a numerical 
code with a specific choice of parameters if the latter is 
far from convergence near the CH.) 

Let us consider first radial geodesies. (The case of 
non-radial geodesies will be treated next.) In the null 
case (to = and p = 0) it is easy to solve the geodesic 
equations (g)-(^). For outgoing geodesies one readily 
finds that the solution is u = const and V = const/F. 
The metric function F can be expanded in [— ln(— V)] , 
despite its non-analyticity in V . To the leading orders in 
[-ln(-V)]- 1 we find that 



F = F (u){l + B[-ln(-F)] 



such that 



^^{i-shW-vx 1 } 



(15) 



(16) 
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The solution for V(r) is then given asymptotically close 
to the CH, to the leading orders in [— ln(— I/)] -1 by 

V(T)=r{l-B[-ln(-r)}- 1 }. (17) 

(Recall that the afhne parameter is given up to a linear 
transformation.) To the leading order in [— In (— t)] 1 
we can thus approximate V(t) w r. Note that al- 
though asymptotically V(t) and V(t) behave like r and 
t = 1, correspondingly, V(t) behaves very differently 
from f = 0. In fact, V(t) diverges as r — > 0. There- 
fore, one can approximate V(r) ~ r only if one is in- 
terested in V(t) itself, or at the most in V{t). This 
approximation is invalid for V(t) or higher derivatives. 
For radial timelike geodesies (m = and p — — 1) one 
uses the finiteness of F iU to find approximately that again 
V ~ const/_F, and consequently one finds the same re- 
sult for V(t). In the case of non-radial geodesies (m = 1) 
one can consider a specific value of the retarded time at 
which the geodesic hits the CH singularity. Than, r, r. u , 
and F can be approximated by their values at the sin- 
gularity. When this is done, the equations for null non- 
radial (to = 1 and p = 0) geodesies become inhomoge- 
neous linear equations. The corresponding homogeneous 
equations are nothing but the equations for the radial 
geodesies, which we already solved. Particular solutions 
for the inhomogeneous equations are easy to generate, 
and one finds that again V(t) is given asymptotically as 
before. The last case is the case of non-radial timelike 
(to = 1 and p = —1) geodesies. In this case the geodesic 
equation becomes (under similar assumptions) an inho- 
mogeneous non-linear equation. Although this equation 
is hard to solve directly, it can be checked that the same 
leading order proportionality of V(t) and r is the solu- 
tion also for this case. We thus find that for all causal 
geodesies, to the leading order in [— ln(— t)] -1 , V(t) is 
proportional to r. 

We next re-express 71 as a function of affine parameter 
(proper time) along radial or non-radial null (timelike) 
geodesies. To the leading order in [— In (— r)] we find 
that in a parallelly-propagated frame 

K(t) oc t- 2 [- In (-t)] ~ 2 "~ 2 . (18) 

A necessary condition for a singularity to be strong in 
the Tipler sense is given by the following theorem |p3| : 
For null (timelike) geodesies, if the singularity is strong 
in the Tipler sense, then for at least one component of 
the Riemann-Christoffel curvature tensor in a parallelly- 
propagated frame, the twice integrated component with 
respect to affine parameter (proper time) does not con- 
verge at the singularity. Specifically, the necessary con- 
dition for the singularity to be Tipler strong is that 

X(t) = /" dr' f cIt"\K(t")\ (19) 



does not converge as r — > 0. It can be readily shown 
that when TZ(r) is integrated twice with respect to r, 
I{t) converges in the limit r — > 0. Consequently, we 
find that a necessary condition for any causal geodesic to 
terminate at a Tipler strong singularity is not satisfied. 
Hence, all causal geodesies terminate at a Tipler weak 
singularity, namely, the singularity is Tipler weak. The 
physical content of this result is that the volume element 
of physical objects remains bounded at the singularity. 
We emphasize that this result is valid everywhere along 
the singularity, in particular at late retarded times where 
the non-linear effects (focusing of the generators of the 
CH and the growth of the blue-shift factors) are crucial. 

III. CONCLUDING REMARKS 

Recently, Ori suggested to define a deformationally- 
strong singularity in the following way. Let A(r) be 
a timelike geodesic with proper time r along it. The 
geodesic A(r) terminates at a deformationally-strong sin- 
gularity at r = if at least one of the following two condi- 
tions holds: (i) A(r) terminates at a Tipler strong singu- 
larity, or (ii) there exists a Jacobi field J(r) for which at 
least one parallelly-propagated tetrad component is un- 
bounded at the limit r — > 0]. This definition is more 
physically- motivated than Tipler 's definition, because it 
classifies a singularity as strong not only when the vol- 
ume element vanishes, but also when the volume element 
diverges to infinity, or there is infinite compression in one 
direction, and infinite stretching in a different direction, 
such that the volume element remains bounded. In fact, 
it can be shown that the failure of the necessary con- 
dition for the singularity to be Tipler strong implies not 
only the boundedness of the volume element, but also the 
boundedness of the Jacobi fields themselves P4I , such 
that objects are not expected to be destroyed also be- 
cause of distortions which preserve the volume element 
or divergence to infinity of the volume element. Conse- 
quently, the singularity we are studying here is weak also 
in the sense of Ori (deformationally weak). (Ori's defini- 
tion does not include null geodesies. However, extended 
physical objects move along timelike geodesies, such that 
this deficiency does not restrict our discussion. It is con- 
ceivable that both Ori's definition and theorem for the 
necessary condition for the singularity to be Ori strong 
can be generalized to all causal geodesies.) 

We note that according to the Krolak classification 
of singularities |^5| this is a strong singularity. Specifi- 
cally, if we integrate over the divergent components of the 
Riemann-Christoffel tensor only once, the integral does 
not converge on the singularity. This means that the ex- 
pansion diverges (negatively) on the singularity (Krolak 
strong), but still the volume element (and the distortion 
in general) remains finite (Tipler and Ori weak). One 
might be worried that even if spacetime were classically 
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extendible beyond the CH, this infinite negative expan- 
sion would inevitably result in unavoidable destruction of 
any extended physical object subsequent to its traversing 
of the CH Of course, any classical extension of ge- 
ometry beyond the CH is not unique. We can, however, 
consider an extension with a continuous (C°) metric and 
a unique C 1 timelike geodesic, and assume that the ob- 
ject follows this geodesic ]27| . Any extension of classical 
geometry beyond the CH (which can be modeled as a 
thin layer wherein the geometry is inherently quantum) 
requires an infinite flux of negative energy traveling along 
the contracting CH. This negative energy flux may then 
act to regularize the expansion, such that the deforma- 
tion rate of physical objects beyond the CH would be 
bounded J27j. (The infinite expansion is likely not to de- 
stroy physical objects up to the CH |^| , in contrast with 
Rcf. p(i| .) Indeed, a simplified two-dimensional quantum 
model shows an infinite ingoing flux of negative energy 
along the CH ||^,^| . More recent semi-classical toy mod- 
els of a quantum field on a mass-inflation background are 
not inconsistent with this picture One should 

not take these quantum results too seriously, however, 
because in these models the semi-classical contributions 
are dominated by the regime where curvature is Planck- 
ian, such that the semi-classical approximation is not ex- 
pected to be valid anymore. Instead, a full quantum 
theory of gravity is of need. Of course, in the absence of 
a valid theory of quantum gravity it is difficult to make 
predictions on the detailed interaction of the thin layer 
of the CH with physical objects, but the evidence we 
currently have do not preclude the possibility of objects 
traversing the CH singularity peacefully. 

I have benefited from useful discussions with Patrick 
Brady and Amos Ori. This work was supported by NSF 
grant AST-9731698. 



APPENDIX 

The independent components of the Ricmann- 
Christoffcl curvature tensor which do not vanish iden- 
tically are: 
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